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Semi-Annual Report 


I. Program Development 

During 1982 work was finished on the computer program development needed 
for our basic non-LTE radiative transfer modeling. As a result it is possible for the 
first time to determine the formation of spectral lines in expanding spherical atmo- 
spheres in a physictMy realistic way, taking into account multilevel atomic 
processes, partial frequency redistribution, and other non-LTE transfer effects tnat 
affect the formation of optically thick lines. This program development work has 
been documented by Avrett and Looser (1983). 

II. Results 

We have begun an exploration of the formation of Mg 11 and Ca I! circumsteliar 
absorption lines in late-type giants and supergiants. This work is based on the 
wave-driven wind theory of Hartmann and MacGregor (1980), as modified by Hart* 
mann, Edwards, and Avrett (1982) and Hartmann, MacGregor, and Avrett (1983) to 
include solution of a realistic energy equation. In this theory, the mass loss rate 
from a star is specified by the chcice of the initial magnetic field strength and 
Alfven wave flux. This choice is guided by observational estimates of stellar mass 
loss rates. In the context of this theory, in which the waves are assumed to dissi- 
pate exponentially with a constant damping length X, the terminai velocity is set by 
the value of X. Although in this schematic treatment there is no theoretical reason 
to choose a specific value for X, it is found that the only way that low terminai 
velocities in agreement with observation can be produced is by setting X ~ 1R*. 
With the wave flux end the damping length fixed, the wave heaNng rate is deter- 
mined. In these winds, the temperature is determined from balancing wave heating 
with radiative cooling. We havs calculated the radiative cooling rate as a function 
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of density and temperature from the results of plane-parallel chromospheric models 
and use these results to approximate the relative cooling in an extended wind 
(preliminary calculations Indicate that this la roughly correct). This permits the run 
of temperature to be calculated along with the density and velocity profiles. 

The most Important prediction of these models Is that a warm zone In the wind 
must exist a result of the wave heetlng. Within this 2 one. the Ca II and Mg II 

atoms can be Ionized to Ca III and Mg III, so that the gas Is transparent In the 
resonance transitions. At large distances, the wind cools adlabaticaily, and the Mg 
III and Ca III can recombine to the lower stage of Ionization. Since the gas Is dif- 
ferentially expanding, this behavior gives rise to "double" clrcumstellar absorption. 
Close to the star, the gas temperature is not high because the high densities le^d 
to efficient radiative cooling. Thus, low-velocity Mg II and Ca II absorption Is 
formed. As the wind expands, the gas heats up and becomes transparent. At 
large distances and high velocities, the wind cools, and the recombination once 
again results In observable clrcumstellar absorption, In this case with a large blue- 
shift. 

The line profiles resulting from this wind behavior are shown in Figs. 1 , 2, and 
3. The atmospheric mod el- the run of density, temperature, and velocity-has been 
computed with a program solving the energy and momentum equations self- 
consistently (Hartmann, Edwards, and Avrett 1982). This model orovides the Input 
to the radiative transfer code. 

The Ionization balance Is solved for hydrogan using a model atom with three 
levels plus continuum. The resulting solution provides the starting polit for calcula- 
tions for a 3-level Ca II model atom and a 6-level Mg II model atom, solving the 
transfer problem in partial redistribution In a moving atmosphere as described 
above. 


The results show the sensitivity of the high- and low- velocity absorption to 
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variations In the mass loss rate and in the maximum wind temperature. Comparison 
with some observed line profiles shows a striking resemblance between theory and 
observation. 

Another feature worthy of notice is that the low-velocity absorption component 
is slightly hlueshifted by ~ -10 to -20 km s’ 1 . A similar shift is observed in many 
stars with high-velocity circumstellar shells (cf. Reimerr 1982). Although there is 
some confusion in the analysis of the low-velocity Mg II and Ca II features due to 
the possibility of interstellar absorption, no such contamination exists for Ha. Figs. 

I and 2 show a similar, small blueshift in the line core; there is a great deal of ob- 
servational evidence for shifts in G and K supergiants (Mallik 1982; Zarro 1982), 
indicating that the flow is beginning in the low chromosphere. 

Note that in the tow mass loss rate case, model 1, the high-velocity absorption 
is much less prominent than in model 2. The high-velocity circumstellar absorption 
in model 1 may be overestimated, because at the radius where recombination is be- 
ginning, the flow times are comparable to the recombination times. This convective 
effect, which has not been included in the statistical equilibrium equations, will tend 
to reduce the Ca It and Mg II absorption. Thus it is clear that minor changes in the 
properties of model 1 could result in no high-velocity absorption at all. The only in- 
dication of mass loss would be the asymmetry observed in Mg II, which is not 
present in Ca II. Stencel and Mullan (1981) have shown that stars exist which 
show no high-velocity absorption, and have Mg II mass-loss asymmetries but nc Ca 

II asymmetry. This is naturally explained by the warm wind model. 

Our results suggest that the warm wind model is correct in outline, which has 
important implications for the physical understanding of late- type stellar winds, it 
may be possible to investigate the required damping length empirics 'ly through cal- 
culations of the type demonstrated here. 


We are continuing an investigation into the wind of a Ori. We showed previ- 
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ously that the warm wind theory predicts Ha and free-free emission in approximate 
agreement with speckle data (Goldberg et a/. 1 982) and radio observations (Newell 
and Hjellmlng 1982). This work Is continuing In an effort to determine wind tem- 
peratures more precisely. Our prior calculations assumed detailed balance in the 
Lyman lines. This assumption Is not adequate for accurate estimation o' the wind 
temperature, and we are recomputing these models. Ultimately we hepe to derive 
the electron temperature and hence the importance of thermal effects in the wind 
dynamics, as well as determine the amount of mechanical energy required to heat 
the wind. 

III. Further Work 

In addition to exploring further the a Or! series of models, we wish to reexam- 
ine the T Tauri series of models (Hartmann, Edwards, and Avrett 1982). These 
models were originally computed with the statistical equilibrium calculated from a 
static model, and we wish to chock the effects of expansion on the source func- 
tions, as well as departures from detailed balance in the Lyman lines. 
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